Recolliding periodic orbits 
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We show that a family of key periodic orbits drive the recollision process in a strong circulary 
polarized laser field. These orbits, coined recolliding periodic orbits, exist for a wide range of 
parameters and their relative influence changes as the laser and atomic parameters are varied. We 
find the necessary conditions for recollision-driven nonsequential double ionization to occur. The 
outlined mechanism is universal in that it applies equally well beyond atoms: The internal structure 
of the target species plays a minor role in the recollision process. 
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Can an ionized electron be driven back to the core 
by an ultrastrong laser pulse? This is a high-stakes is- 
sue in attosecond physics since the returning electron, 
by carrying back the energy it has absorbed from the 
laser, can act as the agent of many key processes in in- 
tense laser physics pQ, including the ultrafast imaging 
of macromolecules [2j and the design of new light sources 
through generation of ultrahigh harmonics [3 . When the 
laser is linearly polarized, the "recollision" (or "rescat- 
tering") |4] [5] model has been immensely successful in 
both interpreting current experiments and devising new 
ones. Called the "keystone of strong-field physics" [6], 
it describes how an ionized electron returns to the core 
in the next half laser cycle to share the energy it has 
acquired from the field with its parent ion. This energy 
may result in the emission of coherent radiation (high 
harmonic generation [3 E]) or, if the energy exchange 
is sufficient, it may alter the core structure, e.g., lead- 
ing to nonsequential double (or even multiple) ionization 
(NSDI) or fragmentation in molecules [9j [10]. For lin- 
ear polarization, an experimental signature of NSDI, and 
thus recollision, is the characteristic "knee" shape in the 
double ionization yield versus laser intensity. The argu- 
ments explaining recollision with linear polarization also 
predict its absence in a circularly polarized (CP) field 
because ionized electrons tend to spiral away from the 
core [TT]. When a knee was found for magnesium with 
a CP field [12] a decade ago, recollision was immediately 
ruled out as a possible explanation. Recently, we rec- 
onciled this surprise with other experimental results in 
CP where no knee was observed [13 , by showing that 
recollision is possible with CP [T4] , 

In this Letter we show that recollisions result from a 
subtle compromise between the action of the strong laser 
field (which leads the electron away in a swirling mo- 
tion) and the Coulomb attraction (which tends to recall 
it). Special families of orbits, which we call recolliding 
periodic orbits (R-PO), turn out to be this compromise: 
They drive the recollision process by repeatedly ioniz- 
ing and returning to the core. Considered together with 
an energy condition to alter the core structure, R-POs 
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FIG. 1. (color online) Left panel: Typical NSDI of Mg as 
modeled by Hamiltonian |2]). The dark (red online) electron 
exhibits a recollision which causes the gray electron to ionize. 
Right panel: Typical one-electron trajectory of Ceo which also 
exhibits a recollision [15]. The gray annulus represents the 
potential well. In both panels, the segment of the trajectory 
which is in bold mimics the shape of key periodic orbits. The 
laser parameters are 780 nm wavelength and an intensity of 
5 x 10 13 W • cm . All trajectories are shown in the rotating 
frame. 

give rise to a few useful rules-of-thumb which can pre- 
dict whether or not a target species (atom, molecule, or 
cluster) will exhibit NSDI, and if so, at which intensities. 

A typical NSDI trajectory in a CP field for a two- 
active-electron atom appears in the left panel of Fig. [I] 
The dark electron (red online) recollides and causes the 
other electron (light gray), previously bound to the nu- 
cleus, to ionize. For comparison, the right panel of Fig.[l] 
displays a recolliding trajectory for Ceo (buckminster- 
fullerene), where the gray annulus represents the poten- 
tial well of the molecule. In both panels, we highlight the 
R-PO signature in bold. We compare these recolliding 
trajectories with the periodic orbits of an already- ionized 
electron in Fig. [2] for each system respectively: The good 
agreement between the shapes of the R-POs of Fig. [2] and 
the sample trajectories of Fig. [I] (bold portions) forms the 
nub of our argument that specific periodic orbits drive 
recollisions. 

The trajectories of Figs, [I] and [2] are represented in 
a frame rotating with the CP field (referred to as the 
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FIG. 2. (color online) R-PO corresponding to the highlighted 
trajectory portions in Fig. [I] The left panel corresponds to 
the hard Coulomb potential, as given by Hamiltonian and 
the right panel corresponds to the one-electron Hamiltonian 
model for Cqq [15] . 

rotating frame). They each consist of an interior loop, 
occurring in the down-field direction of the laser, which 
passes close to the nucleus (where the Coulomb attrac- 
tion dominates the dynamics) and a farther reaching ex- 
terior loop which encloses the nucleus (where the laser 
field predominates). Comparing the R-PO in Fig. [2] for 
Hamiltonian ([TJ (left panel) and the same R-PO for Ceo 
(right panel), we observe striking similarities, despite the 
strong differences in the potentials. The only common 
feature between the two potentials is the Coulomb tail 
far from the core (— 1/r for r ^> 1) which results in a one 
electron model which we examine next. 

In the rotating frame the Hamiltonian modeling a one- 
electron dynamics reads 

p 2 + p 2 1 

£ = - o v / o + E ° x u ( x Py yPx) ' ( 1 ) 

z yr + y z 

where the energy JC is referred to as the Jacobi constant 
for its link to celestial mechanics [16]. The variables (x, y) 
and (p x ,Py) are the position and canonically conjugate 
momentum of the electron. Eq is the amplitude of the 
field and uj its frequency. For Ceo we use a continuous 
approximation of the potential given in Ref. [15] . The 
phase space of Hamiltonian is unbounded, but not 
all electrons can leave the core region and ionize. In 
particular if the Jacobi constant of the electron is smaller 
than the one corresponding to a specific Stark saddle 
point [T7HT9] . the electron is stuck in the core region 
with no possibility for ionization. In Fig. [3j we display 
the limits of the domain accessible to the electron (in 
configuration space), where the saddle point is indicated 
by a sphere. 

As illustrated by Figs. [I] and [2] recollisions with CP are 
organized by certain types of periodic orbits. Only those 
periodic orbits which contain segments both close to and 
far away from the core play a role in recollisions and these 
are R-POs. The topologically simplest R-PO consists of 
an off-centered circle and we found such an orbit for high 
values of the Jacobi constant (top orbit in Fig.[3|. When 




FIG. 3. (color online) Limits of the accessible domain for 
an electron modeled with Hamiltonian |l]). We display three 
R-POs of the family O n - The saddle point is marked by a 
sphere. The laser parameters are the same as in Fig. ^ 

followed by continuity as the Jacobi value is decreased, 
this orbit bifurcates into a family of periodic orbits which 
consist of an interior loop, in the down- field direction of 
the laser, and one or several exterior loops, leading to 
an excursion far from the core (like the highlighted re- 
gions of Fig. [I]). We label these R-POs (9 n , where n 
corresponds to the number of loops in the periodic or- 
bit (see Fig. |3|. The interior loop is responsible for the 
exchange of energy between the already-ionized electron 
and the core (or second) electron while the exterior loop 
is where the ionized electron gains energy from the laser 
field. Not all R-POs are equally important for recollision: 
Relevant orbits have a period which is much smaller than 
the pulse duration (so as to influence the motion of the 
already-ionized electron), and are weakly hyperbolic, so 
that an electron can stay close to them long enough to 
imitate their dynamics. We found a handful of such R- 
POs and their influence waxes and wanes with the choice 
of Jacobi value and intensity. For instance, the influence 
of O4 (middle curve in Fig. |3| can be seen on Fig. 2 of 
Ref. P3]. 

The main outcome of the recollision is the modification 
of the core structure, leading to, e.g., NSDI. Recollisions 
exhibited by the two-electron Hamiltonian follow the or- 
ganizing structures of the one-electron Hamiltonian (as 
seen in Fig. [T]) and indeed it will be shown that these 
structures dictate the properties of the NSDI channel. 
The two-electron dynamics can be expressed by [20] 

^ = iPii 2 + | P2 | 2 +y(|ri|) + y(|r2|)+ 1 

2 v /| ri _ r2 | 2 + 6 2 

Here i*i ? 2 and pi 5 2 are the canonically conjugate positions 
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FIG. 4. (color online) Curves: Double ionization yield for 
Mg (left panel) and Xe (right panel), obtained from Hamilto- 
nian |2|. In both panels the gray shaded background corre- 
sponds to the probability that the three criteria for NSDI are 
satisfied. The laser wavelength is 780 nm. 



and momenta of the two electrons in the lab frame. The 
potential is chosen as V(r) — —2/\Jr 2 + a 2 , where a is 
the electron-core softening parameter [2TJ [22] which is 
adjusted to model the various atoms under investigation 
(for Mg, we set a = 3, for Ne, a = 1, for Ar, a = 1.5, 
and for Xe, a = 1.8). The pulse envelope is given by 
/. For the numerical results presented here, we consider 
a laser envelope with two laser cycle ramp-up and six 
cycle plateau. We find qualitatively similar results with 
other laser envelopes. Initial conditions of the respec- 
tive atoms are taken from a microcanonical distribution 
on the ground state energy (defined as the sum of the 
two first ionization potentials. In this Letter we restrict 
Hamiltonian Q to two spatial dimensions (fully three 
dimensional calculations follow the same organizational 
structure). In Fig. [4] we show the probability of double 
ionization for Mg (left panel) and Xe (right panel) , where 
we use a distance criterion for ionization. Both atoms ex- 
hibit NSDI, which manifests itself in the knee enhance- 
ment. It should be noted that for Ar and Ne, such knees 
are absent in the double ionization probability curve ver- 
sus intensity. In all cases, trajectory inspection shows 
that NSDI corresponds to recolliding trajectories. 

In order to interpret the presence of a knee enhance- 
ment in the double ionization probability versus laser in- 
tensity, we consider the recollision mechanism driven by 
R-POs. Given that R-POs are well defined for a constant 
laser envelope, we consider / = 1 and address recollisions 
happening during the plateau for the two-electron calcu- 
lations. It should be noted that the analysis extends to 
events happening during the ramp-up, by considering a 
lower effective laser intensity. In order to be influenced 
by a R-PO, a trajectory should get close to it and there- 
fore have a Jacobi value compatible with those of the 



R-PO family. The analysis of the recollision dynamics 
shows that the existence of an overlap between the Ja- 
cobi values for the pre-ionized electron and the domain 
of existence of the R-POs accurately predicts the exis- 
tence of recollisions. From this observation, we derive a 
simplified predictive model for the existence of recollision 
and NSDI for a given atom. Given the similarity between 
the R-POs with different models, we use Hamiltonian ([I]) 
to determine the domain of existence of R-POs and re- 
strict the analysis to 02, irrespective of the atom. Due 
to the turn-on of the field the energy gained by each elec- 
tron is E$x^ where is the x-coordinate of the i th 
electron at the beginning of the plateau. The generated 
distribution of Jacobi values corresponds to the colored 
areas in Fig. [5] The gray strip shows the domain of exis- 
tence of the R-PO 02 in the parameter space (/, JC) for 
Hamiltonian ([!]). Following the previous discussion, we 
estimate the probability of recollision as the proportion of 
pre-ionized electrons with the Jacobi values compatible 
with the existence of the R-PO 02- Visually, it corre- 
sponds to the overlap between the gray and respective 
colored regions in Fig. [5] (hatched area). There is signif- 
icant overlap between the surfaces for Mg and Xe and 
the domain of existence of 02, therefore recollisions are 
expected for both atoms (at intensities where there is 
overlap). In contrast, Ne and Ar do not show any over- 
lap so recollisions are not expected for these atoms in the 
near infrared regime. These predictions are confirmed by 
two-electron simulations of Hamiltonian Q for the cor- 
responding atomic models. Other R-POs, e.g., 03 and 
04, can be included in the analysis. However, the results 
are very robust and do not change quantitatively because 
of the strong overlap in the domains of existence of the 
individual n in the parameter space (/, JC). 

Not every recollision leads to NSDI. For that to hap- 
pen, the returning electron needs to bring sufficient en- 
ergy to the core region in order to ionize the bound elec- 
tron while remaining ionized itself. Based on this sim- 
ple picture, we further refine the recollision criterion to 
predict the range of intensities where NSDI, and thus 
the knee in double ionization versus laser intensity, is 
expected. For a given recolliding trajectory, we define 
/Ci (resp. /C2) and JC[ (resp. JC' 2 ) as the Jacobi values 
of the pre-ionized (resp. core) electron before and after 
the recollision and we define rj 2 as the energy exchange. 
Since the core electron is not screened by the pre-ionized 
electron, an effective charge of —2 is used to compute 
its Jacobi value from Hamiltonian ([!]). Assuming elas- 
tic recollisions leading to double ionization, the following 
inequalities hold: 

JC[ — fc\ — rj 2 > /C*, (3a) 
/C 2 = /C 2 +77 2 >/C*, (3b) 

where /C* is the energy of the Stark saddle. Com- 
bining the two equations, we arrive at the condition 



4 




Intensity (W cm 2 ) 



FIG. 5. (color online) Gray surface: Domain of existence of 
the R-PO G 2 for Hamiltonian The gray shading corre- 
sponds to the Greene's residue [23], an indication of the sta- 
bility of C>2- Colored surfaces: Distributions of Jacobi values 
of the pre- ionized electron for Mg (red) and Ne (green). Inset: 
Same as the main panel with Xe (magenta) and Ar (blue). In 
both panels, the hatched areas correspond to the overlap of 
the Jacobi distributions with the domain of existence of Oi- 
The laser wavelength is 780 nm. 

/Ci + /C2 > 2/C* which ensures that the returning elec- 
tron is sufficiently energetic to trigger NSDI. In summary 
there are three necessary conditions for recollision-driven 
NSDI to occur: 

(i) the existence of an R-PO (in a single active electron 
approximation) for Hamiltonian ([I]), 

(ii) an overlap between the distribution of Jacobi val- 
ues of the pre-ionized electron and the domain of 
existence of this R-PO, 

(iii) the pre-ionized electron brings in a sufficient 
amount of energy to free the second electron. 

The first two conditions ensure that recollisions are pos- 
sible, while the third criterion ensures that recollisions 
could lead to NSDI. Varying the intensity and estimating 
the probability of these conditions provides the approx- 
imate intensity range where NSDI is possible. Referring 
back to Fig. |4] we compare these probabilities (gray sur- 
face) for Mg (left panel) and Xe (right panel) with their 
respective double ionization curves, as given by Hamilto- 
nian ([2]). We see that the intensity range at which NSDI 
occurs, i.e., the location of the knee, is predicted well 
by the three conditions. Numerically, Xe shows NSDI at 
intensities approximately one order of magnitude larger 
than in Mg [24] (see also Fig. [4|, in agreement with the 
probability to satisfy the three conditions. These predic- 
tions also agree with experimental findings [9] [T2] , Fi- 
nally, since the double ionization curves of Fig. [4] are 



computed with a ramp-up, the electron experiences an 
effectively lower intensity during the ramp- up. This re- 
sults in double ionization curves which exhibit a cut-off 
intensity for the NSDI channel which is slightly higher 
than what is predicted from the three criteria. 

Our discussion on recollision-driven events boils down 
to a few rules-of-thumb which apply to systems more 
complex than atoms. Molecular recollisions with CP are 
usually attributed to the spatial extent of the system, 
where the pre-ionized electron recollides at a different 
atomic center than the one it originates from [10]. In 
contrast, we have argued above that the excursion of the 
electron is much larger than the size of the molecule (see 
Fig. [T]), and therefore the possibility of recollision in a 
CP field for molecules is not due to their spatial extent 
but mainly to how easily the first electron - the energy 
carrier - can be pre-ionized by the laser. For example, 
Ref. [10 reports a knee for NO and none for N2 using a 
near-infrared CP field. The first ionization potentials of 
N2 and Ar are close, whereas NO resembles Mg. Since 
there are no recollisions for Ar at this wavelength (see 
inset of Fig. [5|, none should be expected for N2. On 
the other hand, recollisions are expected for NO since 
they are plentiful in Mg. Furthermore, since the second 
electron is more tightly bound in NO than it is in Mg, 
these recollisions need to bring back more energy for dou- 
ble ionization and hence experimentally NSDI is seen at 
higher intensity for NO than for Mg [10] [12] . 
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